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ABSTRACT 

The transiting hot Jupiter XO-2b is an ideal target for multi-object photometry and spectroscopy as 
it has a relatively bright (U-mag = 11.25) K0V host star (XO-2N) and a large planet-to-star contrast 
ratio (Rp/R s w 0.015). It also has a nearby (31.21”) binary stellar companion (XO-2S) of nearly 
the same brightness (V-mag = 11.20) and spectral type (G9V), allowing for the characterization 
and removal of shared systematic errors (e.g., airmass brightness variations). We have therefore 
conducted a multiyear (2012-2015) study of XO-2b with the University of Arizona’s 61” (1.55 m) 

Kuiper Telescope and Mont4k CCD in the Bessel U and Harris B photometric passbands to measure 
its Rayleigh scattering slope to place upper limits on the pressure-dependent radius at, e.g., 10 bar. 

Such measurements are needed to constrain its derived molecular abundances from primary transit 
observations. We have also been monitoring XO-2N since the 2013-2014 winter season with Tennessee 
State University’s Celestron-14 (0.36 m) automated imaging telescope to investigate stellar variability, 
which could affect XO-2b’s transit depth. Our observations indicate that XO-2N is variable, potentially 
due to cool star spots, with a peak-to-peak amplitude of 0.0049 ± 0.0007 R-mag and a period of 
29.89±0.16 days for the 2013-2014 observing season and a peak-to-peak amplitude of 0.0035±0.0007 R- 
mag and 27.34 ± 0.21 day period for the 2014-2015 observing season. Because of the likely influence 
of XO-2N’s variability on the derivation of XO-2b’s transit depth, we cannot bin multiple nights of 
data to decrease our uncertainties, preventing us from constraining its gas abundances. This study 
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demonstrates that long-term monitoring programs of exoplanet host stars are crucial for understanding 
host star variability. 

Subject headings: planets and satellites: individual (XO-2b) — methods:analytical — atmospheres — 
radiative transfer - planets and satellites: general 


1. INTRODUCTION 

The exoplanet XO-2b is arguably typical of the more 
extensively measured transiting exoplanets, consider¬ 
ing its basic characteristics and available observations. 
It is a Jupiter-sized body that orbits a KOV star at 
a dist ance of 0.0369±0 .002 AU and a period of 2.6 
days dBurke et al.1 1200711. The re latively bright {V- 
mag == 11.25: iBenavides et al.ll2010f) host star, XO-2N, 
has a southern binary companion, XO-2S, that resides 
too far away (460 0 AU) to affect the planetary orbit 
(1 Burke et al.l [2007). Both stars have relativ ely high 
metallicities and C/O ratios (|Teske et al.lI2013h . as well 
as the same stellar type and brightness. Thus, XO-2S 
provides an ideal photometric reference star for transit 
observations of the exoplanet and host star. 

Ground-based observations of XO-2b at optical wave- 
lcngth s reveal the presence of both sodium and potas¬ 
sium (ISing et all 1201 lal I2012H . Broad-band photome¬ 
try measures the planetary radius at which its atmo- 
sphere becomes optica lly thick in the Sloan z band 
dFernandez et al.l l2009tl a nd in the broad 0.4-0.7 /im 
band (jBurke et al.l l2007f) . A large, multi-platform, 
ground-based study of the XO-2 system finds evidence 
for variability in both XO-2N and XO-2S and records 
fourteen R-band phot ometry primary transits of XO-2b 
(lDamasso et al.!l2015H . Secondary eclipse Spitzer /IRAC 
observat ions suggest that XO- 2b has a weak thermal in¬ 
version ([Machalck et al.l I2009T) . The Hubble Space Tele¬ 
scope (HST) with NICMOS obtained a XO-2b transmis¬ 
sion spectrum at near-IR wavelengths of 1.2-1.8 /im, po¬ 
tentially ind icating the presence of water vapor in its 
atmosphere dCrouzet et al.ll2012lf . 

Molecular abundances are difficult to determine from 
near-IR primary transit measurements because the de¬ 
rived mixing ratios depend on t he radius as s umed 
as a function of pressure le vel dTinetti et al.l 120101 
iBenneke fc Seageill2012l 120131) . iGriffithl (|20l4 ) finds that 
the extinction coefficient depends exponentially on the 
pressure-dependent radius, indicating that high precision 
measurements of exoplanetary radii are essential for the 
interpretation of transit spectra. Note that the derived 
mixing ratios are not sensitive, separately, to uncertain¬ 
ties in the host star’s radius, because primary transits 
meas ure the ratio of the planet-to-star radius dGri ffithl 
120141) . Thus the derived planet’s radius scales to the as¬ 
sumed stellar radius, which we take to be 0.964 Rs U n for 
the host star XO-2N. Gas abundances determined from 
primary transit data are also sensitive to the tempera- 
ture profile of the planet, but to a much lesser extent 
(IGriffithll20l1 . 

Here we present repeated ground-based primary tran¬ 
sit measurements of XO-2b’s atmospheric transmission 
with the University of Arizona’s 61” (1.55 m) Kuiper 
Telescope in the Harris B (330-550 nrn) and Bessell U 
(303-417 nm) photometric bands, where the spectrum is 
devoid of molecular features. If the exoplanet is cloud¬ 
less, the opacities within these bands are established by 


H 2 Rayleigh scattering, and thus the atmospheric den¬ 
sity structure and mean molecular weight, rather than 
atomic and molecular features. The mean molecular 
weight is already well known throug h XO-2b’s density 
dNarita et al.ll201 ll I Burke et al.ll2007li . which indicates a 
H 2 -He based atmosphere. Constraints on XO-2b’s ther¬ 
mal profile are provided by prev ious Spitzer /IRAC ob¬ 
servations (iMachalek et all 12009 ). If clouds are present, 
an upper limit to the radius is obtained, as further dis¬ 
cussed below. 

Our study also includes a nightly monitoring program 
of the host star XO-2N for variability with Tennessee 
State University’s C14 (0.36 m) Automatic P hotoelectri c 
Telescope (AIT) at Fairborn Observatory (iHenrvi 4 999 1. 
These R-band photometric measurements suggest that 
XO-2N is indeed variable, likely the result of star spot 
activity, and cou l d influence XO-2b ’ s transit depth (e.g . 
Pont et ahl 120081 lAgol et all 20101 iCarter et all 120111: 


Desert et al. 20111 ISing et al.l 201 lal iMcCullough et all 

2014 lOshagh et al.l 12014 lDamasso et, al.ll2015fl . How- 
ever, it is unclear how much this variability alters XO- 
2b’s signal as all our Kuiper/Mont4k measurements agree 
to ~1<7. Regardless, XO-2N’s variability prevents us 
from binning the multiple years of Kuiper/Mont4k data 
in order to achieve a higher precision measurement of 
XO-2b’s planet-to-star radius ratio. As a result, we can¬ 
not strongly constrain its molecular abundances. 


2. OBSERVATIONS AND DATA REDUCTION 

Photometric measurements were conducted at the Uni¬ 
versity of Arizona’s 61” (1.55 m) Kuiper Telescope, 
equipped with the Mont4k CCD, a 4096x4096 pixel sen¬ 
sor with a 9.7’x9.7’ field of view. Bessel U-band mea¬ 
surements alone were obtained on 2012 January 5, 2013 
January 21, and 2013 February 24 (UT). Harris B-band 
measurements alone were obtained on 2012 October 29 
and 2014 January 30 (UT). Simultaneous U and B-band 
measurements were taken on 2012 December 10, 2014 
February 12, 2015 January 18, and 2015 February 8 
(UT). Typical seeing was ~1.5”. At the telescope, we 
binned the pixels by 3 x 3 to shorten the readout time 
and achieve a plate scale of 0.43”/pixel. To ensure accu¬ 
rate time-keeping, an on-board clock was automatically 
synchronized with GPS every few seconds throughout the 
observational period. 

To help characterize the out-of-transit baseline, each 
set of observations began ~l-2 hours prior to transit 
ingress and ended ~l-2 hours after egress. The flux 
from XO-2N was measured simultaneously with the flux 
from 8 additional stars on the CCD. Each image is bias- 
subtracted and flat-fielded with 10 flats and bias frames , 
using standard IRAF reduction procedures (lTodvifl993 1. 
Although 8 reference stars were imaged, only the bi¬ 
nary companion, XO-2S, is used as a comparison star 
to map out systematic errors (e.g., airmass), because 
it is nearby (31 . 21”) and of similar stellar type (G9V; 
Benavides et al l [2010f) and brightness (U-mag= 11.20; 
Hog et al.ll2000D to the host star. In addition, we distin- 
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guish a short periodicity in one of the stars within the 
held of view (RA: 7:48:05.44; Dec: +50:15:56.1), thereby 
ruling it out as a reference star. 

To extract the time-varying flux of the target and 
comparison stars, we use the Exop la net Data Reduction 
Pipeline lExoDRPL. iPearson et ahll2014f ). Aperture ra¬ 
dius sizes ranging from 7 to 14 pixels (0.98” to 1.96”) 
with steps of 0.1 pixel (0.014”) are explored to incor¬ 
porate sufficient flux from each star and ensure that no 
contaminating light from nearby stars is present in the 
radial profile. A constant sky annulus with an inner and 
outer radius of 16 and 20 pixels (2.24” and 2.8”), re¬ 
spectively, both larger than the target aperture, prevent 
the inclusion of any background star light. The aperture 
radii that produce the lowest scatter in the raw XO-2N 
light curves were used to create the final light curves. 
This selection process minimizes both the point-to-point 
scatter as well as airmass effects. For each night, the light 
curves of XO-2N are then divided by the light curves of 
XO-2S to remove shared systematic errors, the largest of 
which is airmass. 

The light curves are derived from t he da ta with the an¬ 
alytic equations of lMandcl fc Agol (200 2!) to generate a 
model transit. A Levenberg- Marquardt ( LM) non-linear 
least squares minimization dPress et all I1992T) v ia the 
Inter active Data Language mpfit function ( Markwardtl 
I2009T) provides an initial local fit of the model light 
curves to our data. During the entire analysis, the time 
of mid-transit (T c ), planet-to-star radius ratio R p /R s , 
and out-of-transit normalized flux a are left as the only 
free parameters. The orbital period, inclination, scaled 
semi-major axis, eccentricity, and quadratic limb dark¬ 
ening coefficients remained fixed for this analysis and are 
listed in Table CD In cases where the reduced chi-squared 
(Xr) °f the data to the initial-fit model is greater than 
unity, we multiply the photometric error bars by yfx* 
to co mpensate for the underestimated observat ional er¬ 
rors dBruntt et al.ll2005 iSouthworth et al.ll2007 'l in order 
to produce a global best-fit solution with a reduced \ 2 
equal to unity. This inflation technique captures some 
of the additional non-photon (red) noise in the data; the 
time-correlated component of the red noise is captured 
later in the analysis with a prayer bead. 

To find a global fit solution, we employ 4 sim u ltane- 
ous Markov Chain Monte Carlo (MCMC; iFordl l200 ~> : : 
chains of 200,000 links (steps) seeded by randomly per¬ 
turbing the best-fit parameters from the LM analysis. 
Each MCMC chain draws from a prior Gaussian distri¬ 
bution determined by the initial LM fit and tuned so that 
each fitted parameter has an acceptanc e ratio of ^0.25 
to maximize chain efficiency (lFordll2005li . After running 
each MCMC chain, we define the burn-in point for each 
chain as where the \ 2 value first falls below the median of 
all the x 2 values in the chain, and re move a ll l inks prior 
to th is point. A Gelman-Rubin test dGelman &; Rubiiil 
IT992l! finds a scale redu ction factor <1.001 for all fitted 
parameters ( Fordl 12005 1. suggesting that all chains con¬ 
verge to the same global solutions. 

Since the MCMC draws from a Gaussian prior distri¬ 
bution and assumes that every measurement is uncorre¬ 
lated, it can underestimate the uncertainties in the de¬ 
rived para meters, particula rly if there is non-Gaussian 
red noise dCarter & Winnl 12009T) . Therefore we use a 


residual permutation, or “prayer bead”, method (e.g. , 
Jenkins et al.ll2002t iSouthwortlj 12008 iBean et all 120081 
Winn et al. I2008li to give a more robust estimate of the 
parameter uncertainties. The prayer bead method circu¬ 
larly permutes the residuals around the best-fit MCMC 
solution, and solves for the transit depth, mid-transit 
time, and out-of-transit normalized flux a with a LM for 
each permutation. In this way n (where n is the number 
of data points) new “simulated” lightcurves are formed 
in order to generate new posterior distributions of the pa¬ 
rameters. To remain conservative, we choose the larger 
of the two uncertainties from either the MCMC or the 
prayer bead as the final error bars for each derived pa- 
rame ter and al low ou r uncertainties to be non-symmetric 
(e.g.. [Tbdorov et al.ll2012lh The final model fits to the 
data are presented in Tableland Figures [THU 

3. EFFECTS OF STELLAR VARIABILITY 

Initial analyses of the 2012 U and B-band data indi¬ 
cated 1-a agreement between the derived R p /R s values. 
This consistency seemed to reinforce previous studies 
that concl uded that XO-2 i s a non-active system. For 
example, I Burke et all (f2007l l found that all of their Ex¬ 
tended T eam XO-2b lig ht c urves agreed to ~0.5%. In 
addition, iSing et all (|2011al) found no evidence for oc¬ 
culted starspots in their GTC/OSIRIS spectrophotomet- 
ric study of XO-2b. However we found that our out-of¬ 
transit baseline (as indicated by the out-of-transit nor¬ 
malized flux a; see Table H varied 0.6% in the U and 
2.4% in the B, suggestin g that the host sta r XO-2N is 
variable. For comparison. iBurke et al.l (120071 ! found that 
XO-2N was 93.76% (0.07 ± 0.008 mag) the brightness 
of XO-2S in the B-band, whereas we find that XO-2N 
varies between 92.258% and 93.839% of the brightness of 
XO-2S. 

While we do not see any conclusive evidence for 
starspot crossings in any of our light curves, unocculted 
starspots can cause a transit depth to deepen, resulting 
in the retrieval of an incorrect R p /R s measurement and 


an incorrect interpretation of the exoplanet’s atmosphere 

(e.g., 

Pont et al. 2008 Agol et al.l 

20101 ICarter et all 

20 Si 

Desert et all 120111 ISing et al.i 

2011bl). To inves- 


tigate the possibility of photometric variability in XO- 
2N and XO-2S, we made nightly observations of both 
components throughout the 2013-2014 and 2014-2015 
observing seasons with the Tennessee State University 
Celestron-14 (C14) Automated Imaging Telescope (AIT) 
at Fa i rborn Observatory in southern Arizona (e.g.. lHenrvl 
119991 lEaton et al1l2003f) . We acquired ~170 and ^100 
nightly photometric measurements in the Cousins R- 
band during the two observing seasons. We compute 
differential magnitudes of XO-2N and XO-2S with re¬ 
spect to the mean brightness of several of the most sta¬ 
ble comparison stars in the CCD field of view and then 
perform separate frequency analyses of XO-2N and XO- 
2S for each observing season (Table 3 and Figs. 5-8). A 
few outliers were removed in each analysis, based on 3- 
<7 filtering. Further details of our C14 data acquisition 
and reduct i on p rocedures can be found in the paper by 
ISing et ahl (120151 ). which describes a similar analysis of 
the exoplanet-host star WASP-31. 

The 2013-2014 observations of the northern compo¬ 
nent, XO-2N, have more scatter than the 2014-2015 
season (Table 3 and Figs. 5 & 6). The frequency 
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Table 1 

Fixed Model Values 


Parameter 


Values 


Reference 


XO-2b’s Orbital Parameters 

Period (days) 

Inclination (°) 
a/Rs 

Eccentricity 


2.61586178 

88.01 

7.986 

0 


^ing^et^aL (2011a) 
Crouzetg 3 tgd~° ^ 26 l 2 ) 
Crouze^tjffi (2012) 
JJrouzet^tjrL (20L2) 


Host Star XO-2N Parameters 

W// 

log g (cgs) 

[Fe/H] 

Limb darkening coefficients (U-band) 
Limb darkening coefficients (B-band) 


5343 K 
4.49 
0.39 

0.98884840, -0.13270563 
0.81461289, 0.026193820 


> Teske_et = al. (2013) 
i Teske_et i _al. (2013) 
Tcskc ct ak (2013) 
EaslMiaij_e1_aL (2013) 
Eastmai^etjiL (2013) 


Please note that thelTe ske et ah ( 2013 ) values used here agree with a more contemporary 
study bv IDamasso et all (20 i op to lcr. 


Table 2 

Derived Parameters 


Date 

(UT) 

Rp /Rs 

Uncertainty 

Mid-transit time 

(T c - 2456000 JD) 

Uncertainty 

a t 

Uncertainty 

& / &photon 

Bessel U-band 

2012 Jan 5 0.1037 

-0.0035 

+0.0041 

-68.2337 

-0.0030 

+0.0016 

0.88703 

-0.00032 

+0.00041 

3.71 

2012 Dec 10 

0.1047 

-0.0024 

+0.0028 

271.8283 

-0.0025 

+0.0014 

0.88254 

-0.00024 

+0.00024 

3.17 

2013 Jan 21 

0.1035 

-0.0019 

+0.0023 

313.6810 

-0.0018 

+0.0024 

0.93235 

-0.00014 

+0.00031 

3.84 

2013 Feb 24 

0.1059 

-0.0015 

+0.0015 

347.6865 

-0.0006 

+0.0006 

0.87921 

-0.00023 

+0.00023 

3.90 

2014 Feb 12 

0.1072 

-0.0020 

+0.0020 

700.8276 

-0.0012 

+0.0017 

0.87379 

-0.00029 

+0.00028 

3.36 

2015 Jan 18 

0.1031 

-0.0035 

+0.0045 

1040.8883 

-0.0022 

+0.0039 

0.87210 

-0.00019 

+0.00033 

2.79 

2015 Feb 8 

0.1087 

-0.0051 

+0.0033 

1061.8169 

-0.0040 

+0.0023 

0.86897 

-0.00021 

+0.00034 

2.98 

Harris B-band 

2012 Oct 29 0.1054 

-0.0015 

+0.0016 

229.9736 

-0.0005 

+0.0010 

0.93839 

-0.00014 

+0.00027 

3.64 

2012 Dec 10 

0.1030 

-0.0021 

+0.0018 

271.8288 

-0.0023 

+0.0008 

0.92409 

-0.00018 

+0.00019 

3.83 

2014 Jan 30 

0.1042 

-0.0033 

+0.0075 

687.7504 

-0.0017 

+0.0031 

0.93024 

-0.00050 

+0.00057 

6.08 

2014 Feb 12 

0.1076 

-0.0038 

+0.0039 

700.8291 

-0.0030 

+0.0028 

0.92258 

-0.00039 

+0.00044 

5.20 

2015 Jan 18 

0.1038 

-0.0041 

+0.0023 

1040.8889 

-0.0012 

+0.0034 

0.93782 

-0.00022 

+0.00019 

6.19 

2015 Feb 8 

0.1030 

-0.0012 

+0.0028 

1061.8164 

-0.0007 

+0.0025 

0.93196 

-0.00010 

+0.00019 

3.71 


t Out-of-transit normalized flux, effectively measures the brightness ratio of the host star XO-2N to the comparison (and binary com¬ 
panion) star XO-2S 


spectrum in the middle panel of Figure 5 reveals low- 
amplitude brightness variability in XO-2N with a period 
of 29.9 ± 0.16 days. A least-squares sine fit to the phase 
curve in the bottom panel of Fig. 5 shows a peak-to- 
peak amplitude of 0.0049 ± 0.0007 mag, or 0.53%, a 7a 
detection. We take this to be the rotation period of the 
northern component made visible by the rotational mod¬ 
ulation of starspots across the face of the star. This pe¬ 
riod is consistent with the range of 29-44 days predicted 
from the stellar radius of 0. 97 R ^ and u s ini = 1.4 ± 0.3 
from the discovery paper of lBurke et al.l (l2007f) . Bright¬ 
ness variation in the northern component is reduced in 
the 2014-2015 observing season, but the data essentially 
confirm the results from the first season. The 2.5 day 
difference in the two periods likely results from spot evo¬ 
lution on the timescale of the rotation period. 

The observations of the southern component, XO-2S, 
have very low scatter, 0.0025 and 0.0023 mag in the first 
and second seasons, respectively (Table 3 and Figs. 7 
& 8). These values are consistent with our measure¬ 
ment precision for a single observation, as determined 
by inter-comparison of the most constant stars in the 
CCD frame. Period analyses of XO-2S reveal no sig¬ 


nificant periodicity in either observing season. Our ob¬ 
served levels of variability of XO-2N and XO-2S are con¬ 
sistent with their activity levels, log 77%^, of —4.9 1T 0.01 
and — 5.02 ± 0.01, respectively, given bv iDamasso et all 
(120151) . Howe ver our period analys es are inconsistent 
with those of fPamass o et al.l ([2015), despite their use 
of a similarly-sized telescope (40 cm) relative to our 
study (36 cm) and our two monitoring programs par¬ 
tially overlapping in the 2013-2014 season. This discrep- 
ancy probab ly arises due to the difference in sampling: 
IDamasso et al.l (I2015D monitored the XO-2 system for 
42 nights in the /-band over the 2013-2014 season while 
our program measured it at a much higher cadence over 
171 nights in the 2013-2014 season and 97 nights in the 
2014-2015 season. 

Assuming XO-2N’s variability is constant over the last 
4 years of observation^, we can plot XO-2b’s transit 
depth in each filter vs. XO-2N’s 77-band variability 

1 Please note that since the coverage of star spots varies with 
time, the amplitudes of the variability observed in the 2013-2014 
and 2014-2014 seasons are not necessarily applicable to all of our 
transit observations which were recorded in the 2011-2015 winter 
seasons (Table 2). 
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Orbital Phase 


2014 February 12 - Bessell U 



2012 December 10 - Bessell U 2015 February 8 - Bessell U 




2013 January 21 - Bessell U 
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Orbital Phase 



2013 February 24 - Bessell U 



Figure 1. Bessell U-band light curves from the 2012-2013 observ¬ 
ing seasons showing the transit of XO-2b with their best-fit models 
indicated with the red line. The lcr error bars include the readout 
noise, the poisson noise, and the flat-fielding error. The residuals 
for each light curve are in the lower panel. 


Figure 2. Bessell U-band light curves from the 2014-2015 observ¬ 
ing seasons showing the transit of XO-2b with their best-fit models 
indicated with the red line. The lcr error bars include the readout 
noise, the poisson noise, and the flat-fielding error. The residuals 
for each light curve are in the lower panel. 

across the 2013-2014 and 2014-2015 seasons (Figs.[9]and 
m- In the U-band, the transit depth has a —0.60 cor¬ 
relation with the host star’s 2013-2014 variability and 
a 0.45 correlation with the host star’s 2014-2015 vari¬ 
ability. In the B-band, the correlation is 0.19 with the 
host star’s 2013-2014 variability and —0.68 with the host 
star’s 2014-2015 variability. These plots suggest a corre¬ 
lation between the measured transit depth and host star 
variability. However since all of the Kuiper/Mont4k tran¬ 
sit depths agree to ~1 <t, we cannot quantify how stellar 
variability is affecting the derived transit depth with any 
statistical significance. Since the variation in the mea¬ 
sured transit depth is potentially being influenced by the 
host star’s activity, the multiple transit depth measure¬ 
ments cannot be binned together to achieve higher pre¬ 
cision. Thus we cannot constrain XO-2b’s transit depth 
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Figure 3. Harris B-band light curves from the 2012-2014 sea¬ 
sons showing the transit of XO-2b, following the same protocol as 
Figure fTI 


Figure 4. Harris B-band light curves from the 2014-2015 sea¬ 
sons showing the transit of XO-2b, following the same protocol as 
Figure 2] 


beyond 11.72% in the U-band and 11.17% in the B-band. 


4. INTERPRETATION 

To assess the effects of disparate values of the mea¬ 
sured radius, we consider several radiative transfer mod¬ 
els of XO-2b’s data. We include Rayleigh scatter¬ 
ing, absorption due to water and pressure-induced hy¬ 
drogen absorption. The atmosphere is assumed to 
be cloudless. Since XO-2b has an equilibrium tem¬ 
perature between HD 189733b and HD 209458b, we 
adopt a thermal profile inter mediate between th ose of 
HD 189733b and HD 209458b llMoses et al.ll20Tl . The 
line par ameters for H 2 O are ca lc ulated from the param ¬ 
eters of iRothman et al.1 (I2009T) : iTashkun et al.l (120031) : 
iBarber et al.l ( 20061) , while hot pressure-induced Iff ab- 
sorption derives from lBorvsowl (|2002l) : IZheng fc Borvsowl 


(I1995D . As shown in Figure QTJ uncertainties in the XO- 
2b’s U and B-band light curve depths lead to strong 
uncertainties in the pressure-dependent radius, which in 
turn leads to poor constraints on the composition, even 
in the event of highly precise near-IR measurement s, cur¬ 
rently possible with Hubble/ WFC3 (lGriffithll2014D . Our 
observations of XO-2b indicate that its host star vari¬ 
ability prevents us from combining data from multiple 
nights in order to reduce our uncertainties. These limita¬ 
tions indicate the great need for measurements on differ¬ 
ent platforms, monitoring programs with one instrument, 
and simultaneous observations at different wavelengths, 
which sample the same effective radius (and can then be 
interpreted together). 

5. CONCLUSIONS 

Our observations indicate that XO-2b’s host star, 
XO-2N, is variable with a peak-to-peak amplitude of 
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Table 3 

SUMMARY OF AIT PHOTOMETRIC OBSERVATIONS XO-2N AND XO-2S 


Star 

Observing 

Season 

Nobs 

Date Range 
(HJD - 2,400,000) 

a 

(mag) 

Period 

(days) 

Full Amplitude 
(mag) 

XO-2N 

2013-2014 

171 

56572-56808 

0.0035 

29.89 ±0.16 

0.0049 ± 0.0007 

XO-2N 

2014-2015 

97 

56945-57154 

0.0027 

27.34 ±0.21 

0.0035 ± 0.0007 

XO-2S 

2013-2014 

174 

56572-56813 

0.0025 



XO-2S 

2014-2015 

105 

56945-57154 

0.0023 






Figure 5. Top: 171 nightly Cousins R-band photometric obser¬ 
vations of XO-2N from the 2013-2014 observing season, acquired 
with the C14 Automated Imaging Telescope at Fairborn Obser¬ 
vatory. Middle: Frequency spectrum of the nightly observations 
finds a best period of 29.89 zb 0.16 days. Bottom: The observations 
phased with the best period give a full amplitude of 0.0049±0.0007 
mag, a 7a result. 

0.0049±0.0007 mag in R and a 29.89±0.16 day period in 
the 2013-2014 season and a peak-to-peak amplitude of 
0.0035 ± 0.0007 R-mag and a period of 27.34 ± 0.21 days 
for the 2014-2015 observing season. This host star 
variability potentially influences XO-2b’s transit depths. 
Due to the stellar variability, our data cannot be binned 
to reduce uncertainties. Our study not only demon¬ 
strates that ground-based monitoring of a host star for 
variability is crucial for transit observations, but also 
that visible and near-IR observations must be taken at 
the same time to insure that the same stellar conditions 
are probed. 

The work by R. T. Zellcm, C. A. Griffith, K. A. Pear¬ 
son, J. D. Turner, and J. K. Teske was supported by the 
NASA Atmospheres Program. 

GWH and MHW thank Lou Boyd for his many years of 
APT support at Fairborn Observatory and acknowledge 
support from NASA, NSF, Tennessee State University, 
and the State of Tennessee through its Centers of Excel¬ 
lence program. 

RTZ thanks Mario Damasso, Nikole K. Lewis, Mah- 
moudreza Oshagh, and Mark R. Swain for their helpful 
discussions. 


Figure 6. Top: 97 nightly Cousins R-band photometric obser¬ 
vations of XO-2N from the 2014-2015 observing season. Middle: 
Frequency spectrum of the nightly observations finds a best pe¬ 
riod of 27.34 d= 0.21 days. Bottom: The observations phased with 
the best period give a full amplitude of 0.0035 d= 0.0007 mag, a 50- 
result. 



Figure 7. Top: The nightly Cousins R-band photometric obser¬ 
vations of XO-2S from the 2013—2014 observing season. Middle: 
Frequency spectrum of these observations finds no significant pe¬ 
riodicity between 2 and 200 days. Bottom: Phase plot shows only 
a spurious period for the southern component in the 2013-2014 
observing season. 





















-0.010 : X0-2S (2014-15) 



Figure 8. Top: The nightly Cousins R-band photometric obser¬ 
vations of XO-2S from the 2014-2015 observing season. Middle: 
Frequency spectrum of the nightly observations finds no significant 
periodicity between 2 and 200 days. Bottom: Phase plot shows 
only a spurious period for the southern component in 2014-2015. 


We thank the anonymous referee for their helpful com- 
merits and suggestions. 
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Figure 9. XO-2b’s Bessel U (top) and Harris B (bottom) pho¬ 
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These plots suggest a correlation between the measured transit 
depth and host star variability in the U-band (—0.60) but not in 
the B-band (0.19). Since XO-2N’s variability is likely affecting the 
derived transit depth in both passbands, we cannot bin the multi¬ 
ple nights of data together to achieve higher precision. 
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Figure 10. XO-2b’s Bessel U (top) and Harris B (bottom) pho¬ 
tometric light curve depths (R p /R s ) 2 recorded by the Kuiper 61” 
telescope from 2012-2015 vs. XO-2N’s R-band variability extrapo¬ 
lated from the AIT monitoring over the 2014-2015 season (Fig.[6|. 
These plots suggest a correlation between the measured transit 
depth and host star variability in both the U-band (0.45) and the 
B-band (—0.68). 
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compared to calculated model spectra. Due to the variance in the 
transit depths measured both with Hubble and the 61”, there is not 
a single radiative transfer model that can interpret all of the data 
simultaneously, preventing us from determining XO-2b’s 10 bar ra¬ 
dius (Rio) and placing strong constraints on its water abundance 
[H 2 0] (e.g., see IGriffithl[2014]) . 
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